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For a wide class of low-energy observables the dominant and the next-to-leading effects 
originate entirely from the gauge boson propagators. These corrections are sensitive to 
heavy virtual particles due to quadratic divergences of the corresponding diagrams and due 
to non-decoupling of the heavy particles within the Higgs mechanism for the generation 
of masses. Only recently have the two-loop sub-leading corrections to the myy — mz 
interdependence of order 0{G']^'m^m\) been calculated in an expansion in terms of the 
large top mass Moreover, the two-loop Higgs boson mass dependence associated with 
the fermion loop to Ar has also become available |0] These results allow to reduce 
the theoretical uncertainties originating from higher order effects of the basic physical 
quantities [Q. To improve the accuracy of the calculations, corrections of order 0{G']pm\) 
must be taken into account, which require in particular the calculation of two-loop self- 
energy diagrams with only gauge bosons and massless particles. Keeping in mind this 
physical application, we present here analytical results for all two-loop on-shell self-energy 
master integrals with one mass. Another application of the given results is the calculation 
of the "naive" part of diagrams arising within the asymptotic expansion when the external 
momentum is on the mass shell of the particle with the large mass 

The aim of this letter is the analytical calculation of a full set of two-loop on-shell 
self-energy master integrals with one mass, Eqs.([l|)-(|^). The diagrams of this type are 
presented on Fig.|l[ It can be shown |P, 0], that the full set of master- integrals consists of 
the following ones: Flllll, FOOlll, FlOlOl, FlOllO, FOllOO, FOOlOl, FlOlOO, 
FOOOOl, Villi, VlOOl, Jill (see Fig.p with all indices (powers of propagators) equal 
to 1, and two integrals of the JOll-type: with indices 111 and 112, respectively (for the 
massless integrals see P, H). The integration-by-parts method f^, allows not only to 
reduce diagrams with higher indices to master-integrals, but also to calculate the master 
integrals themselves 0. The main idea for the latter step is the following: a diagram with 
the index (« + l) on a massive line can be represented as the derivative of the same diagram 
with index i on that line w.r.t. its mass squared. This representation is used in certain 
equations obtained by the integration by part method, yielding a differential equation (in 
general partial differential equations) for the initial diagram (with the index i) which can 
be solved analytically. This approach [|T3l , Q - to construct the (differential) relations 
between diagrams - is called the 'Differential Equation Method'. For some particular 
cases (where the number of the masses is small) these partial differential equations can 
be rewritten as ordinary differential equations 0. 

To obtain the finite part of two-loop physical results one needs to know the F-type 
integrals up to the finite part, the V- and J-type integrals up to the e-part, and the 
one- loop integrals up to the e^-part. In the last years various methods have been sug- 
gested and many analytical results for one-, two- and three-loop self-energy diagrams with 
different values of masses and external momenta have been obtained. In particular, the 
e-part of one-loop integrals is given in Ref.JlBl, the finite parts of two- loop self-energy 
diagrams with different internal masses and arbitrary external momentum are given in 
Refs.[|l2|,|I^, ll^-ll^, the e-part of two- loop tadpole and sunset-type integrals at thresh- 
old (pseudothreshold) are calculated in Refs.p3[-P7|. Very efficient and elegant methods 



for numerical calculation of Feynman diagrams have been presented in Refs.||28||,|p9||. In 



Recently the accuracy of these two results has been tested in g. 

'"^One of the first calculations of this type for the Standard Model and QED was performed in Refs.||ll|. 

HI- ..... . rr. 

^ An example where no reduction to an ordinary differential equation was made, may be found in |lq| . 
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many of these cases the authors have obtained one- dimensional integrals as well for ar- 
bitrary masses and momentum. We do not know, however, of attempts to reduce these 
integrals to generalized polylogarithms or other analytic expressions. 

We consider here the following master-integrals in Euclidean space-time with dimen- 
sion N = 4:-2e: 



0'NS{IJ}{i,j,m) 
J{IJIC}{i,j,k,m) 



J d^kid^k2P^^{ki,Xm)P^'\ki-k2,J 
X P^^\k2-p,lCm) 



m) 



y{XJICC}{i,j,k,l,m) = R-^ j d^kid^k2P^^{k2-p,Xm) 



F{ABXJ}C}{a,b,i,j,k,m) 



xP^^\ki - k2,Jm)P^^\ki,lCm) P^^\k2,Cm] 

m^K-^ I d^kid^k2P^^\ki,Am) 

xP^^\k2,Bm)P^^{ki-p,Xm) 
xP^'\k2-p,Jm) p(^)(A;i - fc2,/Cm) 



where 



K 



r(i + g) 



P^^\k,m) = 



{k^ + m^y 



the normalization factor l/(27r)^ for each loop is assumed, and A, B,X, J',)C = 0, 1. 

The finite part of most of the F-type master-integrals can be obtained from results of 
Ref. in the limit z ^ 1. Nevertheless, to take this limit is not an easy task in general. 



We have used this approach to obtain the results for FlOllO and FOOIOI. The finite 
part of FlOlOl is given in Refs.|3> FOOOOl and FlOlOO are presented in Ref.[|g. 
These integrals we have verified by taking the above mentioned limit z — >■ 1 of results in 
Ref . ||24|| . Film is calculated in Ref.|2^ and we have checked it numerically. Instead of 
the usually taken FOllOl integral [0, we consider Jill as master integral. We recall 
the results of all master integrals for completeness. The last unknown master integral then 
remains to be FOOlll. This appears to be one of the most complicated integrals and 
below we present details of its calculation. 

The finite part of the integrals of V- and J-type can be found in Refs. [^-[^. The 
calculation of the e e (e^) parts of master integrals of this type have been performed by 



the differential equation method [|T^, ^ 



The results for F-type master-integrals are follows: 



F{^i3JJ/C}(l, l,l,l,l,m; 
and the coefficients {aj} are given in Table I: 



aiC(3) + a2^S2 + a3inCi2) + 0{e) 
V3 
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TABLE I 




Finn 


FOOlll 


FlOlOl 


FlOllO 


FOllOO 


FOOlOl 


FlOlOO 


FOOOOl 


ai 


-1 





-4 


-1 





-3 


-2 


-3 


(32 


9 
2 


9 


27 

2 


9 


27 

2 


27 

2 


9 





^3 








1 
3 





1 


1 


'2 
3 


1 



where |Tl|, IT7 



^2 = -^Ch (-] = 0.260434137632 • • • . 
9^3 \SJ 

Here we used the — ie prescription. The results for the remaining master integrals are 
the following ones: 



+ _ In 3 + - + ha2) - fegClS) - 12-^ - 63^2 + b,a2) ln3 + ^^82^ 

2tt 



63^ , Vr , 1 TT 2 &5 TT , 2lLs3( — 

+— 52ln3 + 4^ n3 ^ n23-— 2 -- ^ 

4 2^3 2 v^^^ ^ 2 



(2) 



where the coefficients {hi} are listen in Table II: 





TABLE 


// 








hi 


h2 


h. 






Villi 


-1 


-6 


9 
2 


4 


9 


VlOOl 


3 


8 


3 
2 





21 



Jlll(l,l,l,m) 



17 59 



65 



2£^ Ae 



16 



1117 
^2~ 



52C(2) + 48C(2) In 2 - 28C(3) \ + 0{e" 



(3) 



1 -4£ 



/I TT 2 

^ + 2^ - 



J011('l,l,2,m) = — —P -I — + 2-^ - -C(2) 

' ' ' ^ 2(l-2e)(l-3£) 73 3^^^ 



TT 



+^18-^ - ^C(2) - 6^ + ^C(3) + 2752I + 0{e^) 



(4) 
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J011(l,l,l,m) 



m 



15e 



1 3 vr 

+ 



2 (l-2£)(l-3£)(2-3e) 2^3 



(45 n 9 TT , 81^1 J ./X 1863^ 867 tt 

m 3 + m 3 - in^ 3 - 21^C(2) - 

8 4 ' 4 V3^^ ^ 2 



(5) 



ONSll(l,l,m) 



l-2e 



TC 



TC 



^ + e<^\n3-9S2 



TT 



3-^C(2) K (^(^'^ 



(6) 



where 31 



-2.14476721256949- 



(to our knowledge Lss has appeared for the first time in the calculation of the e-part 



of two-loop tadpole integrals in Ref.p5| 



The expansion of Jlll(l,l,l,m) up to e 



terms is given in Ref. [30 



The above results were checked numerically. Fade approximants were calculated from 

The Taylor coefficients were 



the small momentum Taylor expansion of the diagrams |32 



obtained by means of the package with the master integrals taken from |T^. In this 
manner high precision numerical results for the on-shell values of most of the diagrams 
were obtained pi]. Further we made use of the idea of Broadhurst ||35[ to apply the 



. Further we made use of the idea of Broadhurst 
FORTRAN program PSLQ |^B[ to express the obtained numerical values in terms of a 
'basis' of irrational numbers, which were predicted by our analytical calculation, i.e. the 
differential equation method. 

Let us point out that the numbers we obtain are related to polylogarithms at the sixth 
root of unity and hence are in the same class of transcendentals obtained by Broadhurst 
[ p5|] in his investigation of three-loop diagrams which correspond to a closure of the two- 
loop topologies considered here. 

To demonstrate our method of calculation we present in what follows some details 
for the FOOlll integral. The basis of our approach is the differential equation method, 
i.e. we write a differential equation for FOOlll, the inhomogeneous term of which is 
expressed in terms of some 'simpler' integrals. Before turning to the final step, we discuss 
first of all these latter integrals. All calculations are performed off shell and only at the 
end we take the on-shell limit. The following extra notations are needed: 
a) tadpoles: 



Ti(«,m) 



(m 1 



(47r)Tr(i) 



^The definition is: Lss (x) = — \n 1 2 sin 1 1 d9 
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T2{i, J, k,m) = d^kid^k2P^'\ki,m)P^^\k2,m)P^^\ki- k2,m), 



b) loops with one massive line: 



L{i,j,m) = J d^k P^'\k, 0)P^^\k-p, 



c) sunset diagrams with two massive lines: 



3{k,i,j,m) = J d^kid^k2P^''\k2-p,0)P^'\ki,m)P^^\ki-k: 



d) two-loop diagrams with four propagators: 



(7) 



Viii,j,k,l,m) = J d^kid^k2P^'\k2-p,m)P^^\ki-k2,0)P^''\ki,m)P^^\k2,0), 
V2(i,j, k,l,m) = J d^kid^k2P^^{k2-p, 0)P'^^\ki- k2,m)P'^''\ki,m)P^^\k2,m), 

e) two-loop diagrams with five propagators: 

F{a,b,i,j,k,m) = J d^kid^k2P^^\ki,0)P^^\k2,0) 

xP^'\ki -p,m)P^^\k2-p,m)P^^\ki - k2,m). 

1. One-loop diagrams with one massive line. 

Because it is ultraviolet finite, the integral with shifted indices, L(l,2,m), is more suitable 
for analytic calculations than the 'usual' master integral with all indices equal to one. 
The recurrence relations 



L(l,l,m) ■ (A^-3) = -2m2L(l,2,m) - (p2 + m2)L(2,l,m), 
L(l,l,m) ■ (A^-3) = Ti(2,m) - + m2)L(l,2,m), 



(8) 



allow to express the needed integrals in terms of the diagram L(l,2,m), which can be 
written as: 



L(l,2,m) 



K fi 



ds- 



1 



p2 7o (1 - s)^(l + as)(i+=) 
— log (1 + a) - £ 



p- 



log2(l + a) + Li2(-a) 



\oi (1 + a) 



+2 log (1 + a)Li2(-a) + 2Si,2(-a) - Li3(-a) 



3 



where a = p^ /m? (a = -1 corresponds to the on-shell case). To obtain this result we used 
the representation (see e.g. [|^, |3j 
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J <fkP^\k,nii)P^^\k 



p, ma J 



r{t + j-N/2) 
(4vr)fr(or(j) 



X 



(9) 



2. r/ie sunset diagram with two massive lines. 

The set of master-integrals of this type consists of two integrals, e.g. J011(l,l,l,m) and 
J011(l,l,2,m). To avoid ultraviolet singularities, however, we choose again integrals with 
shifted indices: J(l, 2, 2, m) and J(l, 1, 3, m). Appropriate recurrence relations connecting 
these integrals need to be derived (see also Refs.[^ 0)- Applying the integration by parts 
relation with the massless line as distinctive one , we get 

J(l, l,2,m) ■ (A^-2) 

= 2 J d^kid^k2P^^\k2 - p, 0)P'^^\ki, m)P^^\ki - k2, m) ■ (fca - k^Y{k2 - pY 

= J d^kid^k2P^^\k2-p,0)P'-'^\ki,m)P^^\ki- k2,m) ■ k^{k2-pY, (10) 

where for the last step the symmetry of the integral is taken into account. Expanding the 
scalar product 2{k2, ^2 — p) = + (j9 — ^2)^ — p"^, we have 

J(l, 1, 2, m)-{N-2) = ^Tj{2, m) - ^ J(l, 2, 2, m) 

+^ J d''k,d''k2P^'\k2-p,0)P^^\ki,m)P^'\k^ - k2,m) ■ kl (11) 

Applying now the integration-by-parts relation with the massive line with index 2 as 
distinctive one, we obtain 



J(l, 1, 2, m){N - 6) = -4m2j(l, 1, 3, m) - 2m2j(l, 2, 2, m) 
- J d^kid^k2P^^\k2 - p, 0)P^^\k,, m)P^^\ki - k2, m) ■ kl 

Combining Eq.(|l^) and Eq.(12) we get the final relation: 



(12) 



J(l, 1, 2, m)2(l - 3e) = T?(2, m) - (/ + 2m2) J(l, 2, 2, m) - 4m2j(l, 1, 3, m). (13) 
The second recurrence relation, 



J(l,l,l,m) (1-25) [1-- 



y - (2 - 5e) rrif 



J(l,l,2,m) 



-m*J(l, 2, 2, m) - m\p' + 2m^) J(l, 1, 3, m). 



(14) 



^In the integration- by-parts procedure one line is considered 'distinctive', i.e. its propagator contains 
nothing but the momentum fc^ which is used for the partial integration (see, e.g. the discussion in 
Ref.H). 
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can be obtained by applying the operator (d/dp^^ to the sunset integral J(l,l,l,m), 
where the replacement /c2 — *■ k2+p is performed in (0). Further use is made of the relation 

(a/ap^) V^)(p,m) =Ai{i + l- N/2)P^'+'^\p,m) -Ai{i + l)m^P^'+'^\p,m), 
and the dimension property of the integral (see e.g. in Ref. 



d d \ 

P~— + ^^T^ I ^{hhk.m) = {N -i-j - k) 3{i,j,k,m). 



dp"^ 



is taken into account. The third recurrence relation, which we need, has the following 
form: 



2m^ 

J(2, 1, 1, m) = J(l, 1, 2, m) - ^J(l, 1, 3, m) 



(15) 



Eq.([T5|) is obtained in a similar manner as Eq.(|T^: the operator [d/dp^j is applied to the 
sunset J(l, 1, l,m) integral with the initial distribution of momenta (|^), and to the same 
diagram with the replacement A;2 ^ ^2 + P- Their difference yields the recurrence relation 
(|T3|). We would like to note that Eq.(|TB|) shows that the sunset diagram J(l,l,3,m) 
should be known up to terms of order e^. 

After tedious calculation using representation (|]), we obtain 



p2 



J(l,2,2,m) = ^(log2y + 2£ 



1 



log^ y - C(2) log y + log (1 - y) log^ y 



+21ogi/(3Li2(-i/) + 2U2{y)) - 3C(3) - 12Li3(-y) - 6U,iy) 

1 /I 



J(l,l,3,m) . 

2(1 — 2e) \e 



log ?/ + £:< 2 logy log 



y 



[i + yni-y) 



-C(2) - 6Li2(-y) - 2Li2(y)| + £'|6Si,2(y') - 8Si,2(y) + 24Si,2(-y) 

,2 



-8Li3(y) - 24Li3(-y) - 11C(3) - 2C(2) log 



y 



U2{y) + 3Li2(-?/) 



+4 

+ ^Wl/ + ^log?/log^^^^^^, 



{l + yf{l-y) 
log ((1 + yf{l - y)) + 2 logy log (1 - y) log (y{l - y 

y \. f y 



log 



'l + yf{l-yy 



_3a + 4 
2a 



log^y - el + 24Li3(-i/) + 6C(3) 



+2 log y(c(2) - 4Li2(y) - 6Li2(-y)) - 2 log' y log {I - y) - \og^ y 



+ 0(5^), 



where = (4 + a)/a and ?/ = (r — l)/(r + 1). Putting a = — 1, we get 
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J011(l,2,2,m) 



C(2)-£C(3) + 0(e^ 



(16) 



and 



J011(l,l,3,m) 



_i_ f 1 _ ^ + ,„3 - i<Hi - His; 



+.4 ^ - 27i^ - 14^C(2) - ^-^ In^ 3 + ^5, In 3 K Oie^) ) . (17) 



I 3 
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Note that the integrals J011(l,l,3,m) and J011(l,2,2,m) are much simpler than the 
master integrals (|) and (^). 
3. The FOOlll diagram. 

The diagram FOOlll demands special consideration. The investigations are similar to 
those done in [^. Applying the integration-by-parts relation three times to one of the 
triangles in F(l, 1, 1, 1, 1, m) (every time with different distinctive line) and using the 
identity 



2F(l,l,2,l,l,m) + F(l,l,l,l,2,m) 
we obtain the differential equation: 



d 



dm' 



■F(l,l,l,l,l,m), 



(A^ - 4){p* + 2p^m^ + 3m*) - 2m^(/ + p^m^ + m 



dm? 



F(l,l,l,l,l,m) = /, (18) 



where the inhomogeneous term / contains only diagrams of Vi, J and L-type: 



/ = (p4 _ ||^L(1, 2, m) + L(2, 1, m))L(l, 1, m) 

-Vi(2,l,l,l,m)-V2(2,l,l,l,m) 
+3m2(p2 + m^)V2(l, 2, 1, 1, m) - m'^{p^ - m^)Yi{l, 2, 1, 1, m). (19) 

Expressing all L's in terms of L(l,2,m) (see (|^)) and using the recurrence relations for V- 
type integrals obtained via integration by parts, we derive a new, simpler representation 
for the inhomogeneous term /: 



{p^ — m^) 
(1 - 25) 



I -As 
1 - 2e 



T?(2,m)- 



p2(3 -8£) +m2(l -4£) 



Ti(2,m)L(l,2,m) 



-2(/ + m^) J(l, 2, 2, m) + 4m^J(l, 1, 3, m) + 2mVi(l, 1, 2, 1, m) 



+Qm 



d 



l + {p'-m^)-— m2V2(l,l,2,l,m) 



(20) 



where for the V-type diagrams the following differential equations hold: 



8 



dm? 



Vi(l,l,2,l,m 



= J(l, 2, 2, m) + 2(/ - m^) 



Ti(3,m)L(l, l,m) - J(l,l,3,m) 
V2(l,l,2,l,m) 



(21) 



p^J(l, 2, 2, m) + 3m2T2(l, 2, 2, m), (22) 



and (see e.g. Eq.(5.5) in Ref. 



T2(l,2,2,m) = — .3^2 + 0(e). 



We integrate Eqs.(|l8D, (0), ( ^2]) with the boundary condition that all diagrams tend to 
zero as m^ oo. Finally for the on-shell case we obtain 



F00111(l, 1, 1, 1, l,m) = 8 / arctan 



x(4 - x) \ 4:-x 



f 2x-l\ , , 
X < arctan ^— — In 1 — x + x 

li-x l6 V3 ; 



This integral is evaluated numerically with a precision of 40 decimals. Then the program 
PSLQ (see above) is applied with 15 basis elements occurring in (|1|) to (^ with the result 
given in Table I. 
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Figure 1: The full set of two- loop self-energies diagrams with one mass. Bold and thin 
lines correspond to the mass and massless propagators, respectively. 
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